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Theory of helimagnons in itinerant quantum systems 
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The nature and effects of the Goldstone mode in the ordered phase of helical or chiral itinerant 
magnets such as MnSi are investigated theoretically. It is shown that the Goldstone mode, or 
helimagnon, is a propagating mode with a highly anisotropic dispersion relation, in analogy to 
the Goldstone mode in chiral liquid crystals. Starting from a microscopic theory, a comprehensive 
effective theory is developed that allows for an explicit description of the helically ordered phase, 
including the helimagnons, for both classical and quantum helimagnets. The directly observable 
dynamical spin susceptibility, which reflects the properties of the helimagnon, is calculated. 

PACS numbers: 75.30.Ds; 75.30.-m; 75.50.-y; 75.25. +z 



I. INTRODUCTION 

Ferromagnetism and antiferromagnetism are the most 
common and well-known examples of long-range mag- 
netic order in solids. The metallic ferromagnets Fe and 
Ni in particular are among the most important and well- 
studied magnetic materials. In the ordered phase, where 
the rotational symmetry in spin space is spontaneously 
broken, one finds soft modes in accord with Goldstone's 
theorem, namely, ferromagnetic magnons. The latter 
are propagating modes with a dispersion relation, or 
frequency-wave vector relation, fl ~ fe 2 in the long- 
wavelength limit. In antiferromagnets, the correspond- 
ing antiferromagnetic magnons have a dispersion rela- 
tion f2 ~ |fe|. In rotationally invariant models that ig- 
nore the spin-orbit coupling of the electronic spin to the 
underlying lattice structure these relations hold to ar- 
bitrarily small frequencies O and wave vectors fe. The 
lattice structure ultimately breaks the rotational sym- 
metry and gives the Goldstone modes a mass. In ferro- 
magnets, the low-energy dispersion relation is also mod- 
ified by the induced magnetic field, which generates a 
domain structure. These are very small effects, however, 
and magnons that are soft for all practical purposes are 
clearly observed, directly via neutron scattering, and in- 
directly via their contribution to, e.g., the specific heat. 1 
These observations illustrate important concepts of sym- 
metries in systems with many degrees of freedom with 
ramifications that go far beyond the realm of solid-state 
physics. 2:3 

In systems where the lattice lacks inversion symme- 
try additional effects occur that are independent of the 
spin rotational symmetry. This is due to terms in the 
action that are invariant under simultaneous rotations of 
real space and M, with M the magnetic order param- 
eter, but break spatial inversion symmetry. Microscop- 
ically, such terms arise from the spin-orbit interaction, 
and their precise functional form depends on the lattice 
structure. One important class of such terms, which is 



realized in the metallic compound MnSi, is of the form 
M ■ (V X M). 4 -' 5 They are known to lead to helical or 
spiral order in the ground state, where the magnetiza- 
tion is ferromagnetically ordered in the planes perpen- 
dicular to some direction q, with a helical modulation 
of wavelength 2n/\q\ along the q axis. 6,7 In MnSi, which 
displays helical order below a temperature T c w 30 K at 
ambient pressure, 2n/\q\ « 180 A. 8 Application of hydro- 
static pressure p suppresses T c , which goes to zero at a 
critical pressure p c ~ 14kbar. 9 

In addition to the helical order, which is well under- 
stood, MnSi shows many strange properties that have 
attracted much attention lately and so far lack explana- 
tions. Arguably the most prominent of these features is 
a pronounced non- Fermi- liquid behavior of the resistivity 
in the disordered phase at low temperatures for p > p c . 10 
In part of the region where non-Fermi-liquid behavior is 
observed, neutron scattering shows 'partial' magnetic or- 
der where helices still exist on intermediate length scales 
but have lost their long-range directional order. 11 Such 
non-Fermi-liquid behavior is not observed in other low- 
temperature magnets. Since the helical order is the only 
obvious feature that sets MnSi apart from these other 
materials it is natural to speculate that there is some 
connection between the helical order and the transport 
anomalies. In this context it is surprising that some ba- 
sic properties and effects of the helically ordered state, 
and in particular of the helical Goldstone mode, which 
we will refer to as a helimagnon in analogy to the fcr- 
romagnons and antiferromagnons mentioned above, are 
not known. The purpose of the present paper is to ad- 
dress this issue. We will identify the helimagnon and 
determine its properties, in particular its dispersion rela- 
tion and damping properties. We also calculate the spin 
susceptibility, which is directly observable and simply re- 
lated to the helimagnon. The effects of this soft mode on 
various other observables will be explored in a separate 
paper. 12 A brief account of some of our results, as well 
as some of their consequences, has been given in Ref. 13. 
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One of our goals is to develop an effective theory for 
itinerant quantum hclimagnets. We will do so by de- 
riving a quantum Ginzburg-Landau theory whose coeffi- 
cients are given in terms of microscopic electronic corre- 
lation functions. Such a theory has two advantages over 
a purely phenomenological treatment based on symmetry 
arguments alone. First, it allows for a semi-quantitative 
analysis, since the coefficients of the Ginzburg-Landau 
theory can be expressed in terms of microscopic param- 
eters. Second, it derives all of the ingredients necessary 
for calculating the thermodynamic and transport proper- 
ties of an itinerant helimagnet in the ordered phase using 
many-body perturbation theory techniques. 1 

The organization of the remainder of this paper is as 
follows. In Sec. II we use an analogy with chiral liquid 
crystals to make an educated guess about the wave vector 
dependence of fluctuations in hclimagnets, and employ 
time-dependent Ginzburg-Landau theory to find the dy- 
namics. In Sec. Ill we derive the static properties from 
a classical Ginzburg-Landau theory. In Sec. IV we start 
with a microscopic quantum mechanical description and 
derive an effective quantum theory for chiral magnets. 
We then show that all of the qualitative results obtained 
from the simple arguments in Sec. II follow from this 
theory, with the additional benefit that parameter values 
can be determined semi-quantitativcly. We conclude in 
Sec. V with a summary and a discussion of our results. 
Some technical details are relegated to three appendices. 



II. SIMPLE PHYSICAL ARGUMENTS, AND 
RESULTS 



Hclimagnets are not the only macroscopic systems that 
display chirality, another examples are cholesteric liquid 
crystals whose director order parameter is arranged in a 
helical pattern analogous to that followed by the mag- 
netization in a helimagnet. 14 There are some important 
differences between magnets and liquid crystals. For in- 
stance, the two orientations of the director order param- 
eter in the latter are equivalent, which necessitates a de- 
scription in terms of a rank-two tensor, rather than a 
vector as in magnets. 15 Also, the chirality in cholesteric 
liquid crystals is a consequence of the chiral properties of 
the constituting molecules, whereas in magnets it is a re- 
sult of interactions between the electrons and the atoms 
of the underlying lattice. However, these differences are 
not expected to be relevant for some basic properties of 
the Goldstone mode that must be present in the helical 
state of either system. 16 We will therefore start by using 
the known hydrodynamic properties of cholesteric liquid 
crystals to motivate a guess of the nature of the Gold- 
stone mode in hclimagnets. In Sec. Ill we will see that 
the results obtained in this way are indeed confirmed by 
an explicit calculation. The arguments employed in this 
section are phenomenological in nature and very general. 
We therefore expect them to apply equally to classical 
hclimagnets and to quantum hclimagnets at T — 0, as is 



the case for analogous arguments for ferromagnetic and 
antifcrromagnetic magnons. 



A. Statics 

Consider a classical magnet with an order parameter 
field M(x) and an action 6,7 



S[M] = J 



dx 



-M\x) + -(WM(x)) 2 



+ \ M{x) ■ (V x M(x)) + \ {M\x)) 2 \ . (2.1) 

This is a classical </) 4 -theory with a chiral term with cou- 
pling constant c. Physically, c is proportional to the spin- 
orbit coupling strength gso- The expectation value of 
M is proportional to the magnetization, and it is easy to 
see that a helical field configuration constitutes a saddle- 
point solution of the action given by Eq. (2.1), 

M sp (x) = m (ei cos q ■ x + e 2 sin q ■ x) (2.2a) 
= m (cosg2:,sin<7z,0) . (2.2b) 

In Eq. (2.2a), e\ and e 2 are two unit vectors that are 
perpendicular to each other and to the pitch vector q. 
The chirality of the dreibein {q, ei, e2} reflects the chi- 
rality of the underlying lattice structure and is encoded 
in the coefficient c in Eq. (2.1), with the sign of c deter- 
mining the handedness of the chiral structure. In Eq. 
(2.2b) we have chosen a coordinate system such that 
{ei,e 2 ,q/q} = {x,y,z}, a choice we will use for all ex- 
plicit calculations. We will further choose, without loss 
of generality, c > 0. The free energy is minimized by 
q = c/2a, and the pitch wave number is thus propor- 
tional to <7so- 

Now consider fluctuations about this saddle point. An 
obvious guess for the soft mode associated with the or- 
dered helical state are phase fluctuations of the form 

M(x) = m (cos(qz + (j)(x)), sin(qz + (j)(x)), 0) 

= M sp (x) + m (j)(x) (— sin qz, cos qz, 0) + 0(<j) 2 ). 

(2.3) 

These phase fluctuations are indeed soft; by substitut- 
ing Eq. (2.3) in Eq. (2.1) one finds an effective action 
<Seff[0] = const./ dx (V0(a;)) . However, this cannot 
be the correct answer, which can be seen as follows. 17 
Consider a simple rotation of the planes containing the 
spins such that their normal changes from (0, 0, q) to 
(cti,a2,q), which corresponds to a phase fluctuation 
<j>(x) = a\ x + cc2 y- This cannot cost any energy, yet 



aj + / for this particular phase fluctu- 



ation. The problem is the dependence of the effective 
action on Vj_</>, where V = (V±,d z ). The soft mode 
must therefore be some generalized phase u(x) with a 
schematic structure 



u(x) ~ 4>(x) + X7±tp(x), 



(2.4) 



3 



where <p(x) represents the z-component of the order pa- 
rameter vector M(x). The lowest order dependence on 
perpendicular gradients allowed by rotational invariance 
is V^u, and the extra term in u proportional to Vj_ will 
ensure that this requirement is fulfilled. The correct ef- 
fective action thus is expected to have the form 



S e s[i 



dx 



c z (d z u{x)Y 



- C_L 



(Viu(*))V 



(2-5) 

where c z and cj_ are elastic constants. The Goldstone 
mode corresponding to helical order must therefore have 
an anisotropic dispersion relation: it will be softer in 
the direction perpendicular to the pitch vector than 
in the longitudinal direction. 18 Separating wave vectors 
k = (kj_,k z ) into transverse and longitudinal compo- 
nents, the longitudinal wave number will scale as the 
transverse wave number squared, k z ~ k\jq. The fac- 
tor 1/q 2 in the transverse term in Eq. (2.5), which serves 
to ensure that the constants c z and c± have the same di- 
mension, is the natural length scale to enter at this point, 
since a nonzero pitch wave number is what is causing the 
anisotropy in the first place. A detailed calculation for 
cholesteric liquid crystals 19 shows that this is indeed the 
correct answer, and we will see in Sec. Ill that the same 
is true for helimagnets. 



B. Dynamics 

In order to determine the dynamics of the soft 
mode in a simple phcnomenological fashion we uti- 
lize the framework of time-dependent Ginzburg-Landau 
theory. 20 Within this formalism, the kinetic equation for 
the time-dependent generalization of the magnetization 
field M reads 



d t M(x,t) = -jM(x,t) 



5S 



SM(x) 



- J dyD(x - y) 



SS 



SM(y) 



M(x,t) 

+ C0M), (2-6) 



M(y,t) 



with 7 a constant. The first term describes the precession 
of a magnetic moment in the field provided by all other 
magnetic moments, D is a differential operator describing 
dissipation that we will specify in Sec. II C, and £ is a 
random Langcvin force with zero mean, (£(x,t)} = 0, 
and a second moment consistent with the fluctuation- 
dissipation theorem. 

Now assume an equilibrium state given by Eq. (2.2b). 
In considering deviations from the equilibrium state we 
must take into account both the generalized phase modes 
at wave vector q, which are soft since they are Goldstone 
modes, and the modes at zero wave vector, which are 
soft due to spin conservation. The latter we denote by 
m(x,t), and for the former we use Eq. (2. 3), 21 

M(x,t) = M sp (x) + m(x,t) 

+ mo u(x, t)(— sinqz, cosqz, 0). (2.7) 



The action for u is the effective action given by Eq. (2.5), 
and the action for m is a renormalized Ginzburg-Landau 
action of which we will need only the Gaussian mass term. 
We thus write 

S[m,u] = y J dx m 2 {x) + S cS [u]. (2.8) 

The mass ro of the zero- wave number mode that appears 
here and in the remainder of this section is different from 
the coefficient r in Eq. (2.1), and we assume r > 0. 22 

We now use the kinetic equation (2.6) to calculate the 
average deviations (m(x,t)} and (u(x,t)} from the equi- 
librium state. For simplicity we suppress both the aver- 
aging brackets and the explicit time dependence in our 
notation, and for the time being we neglect the dissipa- 
tive term. With summation over repeated indices im- 
plied, Eq. (2.6) yields 



d t m 3 (x) = -jesij M* p (x) J dy 

SS 



SS 5u{y) 



-7m / dy 



Su(y) 5Mj(x) 
du(y) 



5u(y) 
— sin qz 



cos qz 



dM v (x) 
du(y) 



dM x (x) 



(2.9a) 



and by using Eq. (2.7) in the identity 



5{x — y) = j dz 



Su(x) SMi(z) 



5Mi(z) Su(y) 



we find 



8u(x) 8u(x) 
m I — sin qz + cos qz 



SM x (y) ' -~*~8M y (y)J ^ V) ~ 

(2.9b) 

Using Eq. (2.9b) in Eq. (2.9a) eliminates the integration, 
and using Eqs. (2.8) and (2.5) we find a relation between 
m.3 and u, 



d t m 3 (x) = -7 



5S C 



- 7 {-c z dl + C± _V\/q 2 ) u( 



Su(x) 

A second relation is obtained from the identity 
dtM^x) 



x 



(2.10) 



/ d y E77~r d My). 



5u(y) 

By applying Eq. (2.6) to the left-hand side, and using 
Eqs. (2.8) and (2.7) we obtain 



/ 



dy 



SMijx) 
Su(y) 



O t u(y)="fr a / dy ^ f _ A m 3 {y) 



Su(y) 



or 



d t u(x) = jr m 3 (x). (2-H) 
Combining Eqs. (2.10, 2.11) we find a wave equation 
8 2 t u{x) = - 7 2 r {-c z d 2 z + c ± V 4 ± /q 2 ) u(x). (2.12) 
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This is the equation of motion for a harmonic oscillator 
with a resonance frequency 



0J {k)= 1 rl /2 ^c z kJ+7 ± k\/q 2 (2.13) 



and a susceptibility 



Xo 



(2.14) 



We thus have a propagating mode, the helimagnon, with 
an anisotropic dispersion relation: for wave vectors par- 
allel to the pitch vector q the dispersion is linear, as in 
an antiferromagnet, while for wave vectors perpendicular 
to q it is quadratic. Fluctuations transverse with respect 
to the pitch vector are thus softer than longitudinal ones. 
The nature of the excitation corresponding to the longi- 
tudinal and transverse helimagnon, respectively, is shown 
in Fig. 1. 

For later reference we note that for determining the 
static properties of the helimagnon it sufficed to discuss 
the phase modes at wave vector q, while the dynamics 
are generated by a coupling between the phase modes 
and the modes at zero wave vector. This observation 
gives an important clue for the correct structure of the 
microscopic theory we will develop in Sec. IV. 



TZ/k z 




CO ~ q k z 



2%/k, 



co~ k, 



FIG. 1: Sketch of a longitudinal (k\\q, left panel) and trans- 
verse (fc _L q, right panel) helimagnon. The solid lines de- 
lineate planes of spins pointing out of (dotted circle) or into 
(crossed circle) the paper plane. 



C. Damping 

In order to investigate the damping of the mode we 
need to take into account the dissipative term in the 
master equation (2.6) which we have neglected so far. 
Usually, in the case of a conserved order parameter, the 
damping operator D in Eq. (2.6) is proportional to a 
gradient squared. 20 In the present case, however, one ex- 
pects an anisotropic differential operator, with different 
prefactors for the longitudinal and transverse parts, re- 
spectively. We will see in Sec. IV that in the particular 



model we will consider the transverse part of the gradient 
squared has a zero prefactor. We thus write 



D(x-y)=T5(x-y)d 2 , 



(2.15) 



where T is a damping coefficient. Going through the 
derivation in the previous subsection again, we see that 
Eq. (2.11) remains unchanged except for gradient cor- 
rections to the right-hand side. Equation (2.10), on the 
other hand, acquires an additional term that is of the 
same order as the existing ones, 

d t m s (x) = -7(-c 2 9 2 2 + c ± Vi/g 2 ) u(a;) 

-r Td 2 m 3 (x). (2.16) 

Together with Eq. (2.11) this leads to an equation of 
motion for u given by 

d 2 u{x) = -l 2 ro (-c z d 2 + c ± V 4 Jq 2 ) u(x) 

-r Td 2 d t u(x). (2.17) 

This corresponds to a damped harmonic oscillator with 
a susceptibility 



1 



X 



ujQ(k) — u 2 — iu)j(k) 



where the damping coefficient is given by 
7 (fe) = r Ffc 2 . 



(2.18a) 



(2.18b) 



Now recall that we are interested in systems where the 
magnetization is caused by itinerant electrons. In a sys- 
tem without any elastic scattering due to impurities, and 
at zero temperature, the coefficient T, which physically 
is related to a generalized viscosity of the electron fluid, 
is itself wave number dependent and diverges for k — > 
as T a 1/| fe|. This leads to a damping coefficient 



7(fe -» 0) oc k 2 J\k\. 



(2.19a) 



We see that j(k) scales as k z (for k z ^> k^/q 2 ), or as 

k^J 2 (for k 2 <C k\_ /q 2 ), while the resonance frequency loq 
scales as k z . If the prefactor of the damping coefficient 
is not too large, we thus have j(k) < uJo(k) for all k 
and the mode is propagating. Any amount of quenched 
disorder will lead to T being finite at zero wave vector, 
and hence to 



j(k -> 0) cx k 2 z 



(2.19b) 



In this case, the mode is propagating at all wave vectors 
irrespective of the prefactor. 



D. Physical Spin Susceptibility 

The physical spin susceptibility Xs, which is directly 
measurable, is given in terms of the order-parameter cor- 
relation function. The transverse (with respect to q) 
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components of Xs are given by the correlations of the 
phase 4> hi Eq. (2.3), and are thus directly proportional 
to the Goldstone mode. In a schematic notation, which 
ignores the fact that </> at zero wave vector corresponds to 
a magnetization fluctuation at wave vector ±q, we thus 
expect 



1 



^o(fc) — LU 2 — iLUj(k) 



(2.20a) 



The longitudinal component will, by Eq. (2.4), carry an 
additional factor of fc 2 , and is thus expected to have the 
structure 



fci 



u>o(k) — lu 2 — ilo 7(fe) 



(2.20b) 



Since to ~ k z ~ k 2 ^ in a scaling sense, we see that the 
transverse susceptibility, x^ ~ 1/w 2 , is softer than the 
longitudinal one, xi ~ 



E. Effects of broken rotational and translational 
invariance 

For the arguments given so far, the rotational symme- 
try of the action 5 [A/], Eq. (2.1), i.e., the invariance un- 
der simultaneous rotations in real space and spin space, 
played a crucial role. Since the underlying lattice struc- 
ture of a real magnet breaks this symmetry, it is worth- 
while to consider the consequences of this effect. 

In a system with a cubic lattice like MnSi, the sim- 
plest term that breaks the rotational invariance is of the 
form 6 - 7 

S cuhic [M] = ^ J dx [{8 x M x {x)f + (d y M y (x)) 2 

+ (d z M z {x)f\ . (2.21) 

Other anisotropic terms with a cubic symmetry (see Ap- 
pendix A for a complete list) have qualitatively the same 
effect. In Eq. (2.21), a\ cx g 2 so , with gso the spin- 
orbit coupling strength (see Sec. II A). On dimensional 
grounds, we thus have a\ = bq 2 a, with b a number, and a 
the coefficient of the gradient squared term in Eq. (2.1). 
6^0 leads to a pinning of the helix pitch vector in (1,1,1) 
or equivalent directions (for b < 0), or in (1,0,0) or equiv- 
alent directions (for b > 0). 6 ' 7 In addition, it invalidates 
the argument in Sec. II A that there cannot be a (V_i_(/>) 2 
term in the effective action. However, the action is still 
translationally invariant, so a constant phase shift can- 
not cost any energy. To Eq. (2.1) we thus need to add a 
term 



ocubic 
'-'off 



. bq 2 a f 
U] = —j 



dx (V±u(x)Y 



(2.22) 



This changes the soft-mode frequency, Eq. (2.13), to 



uj (k) 



1/2 
l r 



k 2 z +baq 2 k 2 L +ctk\/q 2 . (2.23) 



We note that, due to the weakness of the spin-orbit cou- 
pling, a q 2 <C 1 , and therefore the breaking of the rota- 
tional symmetry is a very small effect. In MnSi, where 
the pitch wave length is on the order of 200A, while a 1 / 2 
is on the order of a few A at most, the presence of the 
k 2 ± term is not observable with the current resolution of 
neutron scattering experiments, and we will ignore this 
term in the remainder of this paper. 

The above considerations make it clear that the Gold- 
stone mode is due to the spontaneous breaking of trans- 
lational invariance, rather than rotational invariance in 
spin space. Consistent with this, there is only one Gold- 
stone mode, as the helical state is still invariant under a 
two-parameter subgroup of the original three-parameter 
translational group. 23 In this sense, the helimagnon is 
more akin to phonons than to ferromagnetic magnons. 
Let us briefly discuss the effect of the ionic lattice on this 
symmetry, as the helix can be pinned by the periodic 
lattice potential and therefore one expects a gap in the 
magnetic excitation spectrum. To estimate the size of 
the gap, we investigate the low-energy theory taking into 
account only slowly varying modes with \k — q\ < q. In 
a periodic lattice, momentum is conserved up to recip- 
rocal lattice vectors Gj. The leading term which breaks 
translational invariance, M k — > M k e tkr " therefore is of 
the form 

S n = ]T V {klhG .M kl M k2 ...M k jlY,ki-G J \ 

ki,...,k n ,j \ i / 

(2.24) 

where V parameterizes the momentum dependent cou- 
pling strength (and we have omitted vector indices). 
Within the low-energy theory, all momenta are of order 
of q. Therefore, umklapp scattering can only take place 
for n > G/q. In the case of MnSi, where G/q ss 40, one 
therefore needs a process proportional to M 40 to create 
a finite gap! It is difficult to estimate the precise size 
of the gap which depends crucially, for example, on the 
commensuration of the helix with the underlying lattice. 
However, the resulting gap will in any case be unobserv- 
ably small as it is exponentially suppressed by the large 
parameter G/q cx 1/gso- 



III. NATURE OF THE GOLDSTONE MODE IN 
CLASSICAL CHIRAL MAGNETS 

One of our goals is to derive from a microscopic theory 
the results one expects based on the simple considera- 
tions in Sec. II. As a first step, we show that the phe- 
nomenological action for classical hclimagnets given by 
Eq. (2.1) does indeed result in the effective clastic the- 
ory given by Eq. (2.5). A derivation from a microscopic 
quantum mechanical Hamiltonian will be given in Sec. 
IV. 

The classical </> 4 -theory with a chiral term represented 
by the action given in Eq. (2.1) can be analyzed in anal- 
ogy to the action for chiral liquid crystals. 19 In the mag- 
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netic case the chiral term with coupling constant c is of 
the form first proposed by Dzyalishinkski 4 and Moriya, 5 
who showed that it is a consequence of the spin-orbit in- 
teraction in crystals that lack spatial inversion symmetry. 



A. Saddle-point solution 

The saddle-point equation, 8S/5Mi(x) = 0, reads 

(r -aV 2 + cVx +uM 2 (x)) M(x) =0, (3.1a) 

and the free energy density in saddle-point approxima- 
tion is given by 



(3.1b) 



with M sp a solution of Eq. (3.1a). 

The helical field configuration given by Eqs. (2.2) with 
an amplitude m\ = — (r+aq 2 — cq)/u is a solution for any 
value of q. The physical value of q is determined from 
the requirement that the free energy must be minimized, 
which yields 



and 



q = c/2a. 



1/2 



(3.2a) 



(3.2b) 



The zero solution, m = 0, is unstable with respect to 
the helical solution for all r < c 2 /4a, and for c ^ 
the ferromagnetic solution q = 0, Wq = —r/u is always 
unstable with respect to the helical one since one can 
always gain energy by making q ^ due to the linear 
momentum dependence of the chiral term M ■ (V x M). 



B. Gaussian fluctuations 

1. Disordered phase 

In the disordered phase, r > c 2 /4a, the Gaussian prop- 
agator is easily found by inverting the quadratic form in 
Eq. (2.1), 



(Mi{k)Mj{p)) = 5 k ,- P 



1 



(r + aq 2 ) 2 — c 2 q 2 



Sij(r + aq 2 ) + e {j iick t - hkj 



r + aq 2 



. (3.3) 



The structure of the prefactor in this expression is consis- 
tent with the conclusion of Sec. Ill A: For r > c 2 /4a, the 
denominator N(q) = (r + aq 2 ) 2 — c 2 q 2 is minimized by 
q = 0, and N(q) has no zeros in this regime. N(q) first 
reaches zero at r = c 2 /4a and q — c/2a, and the disor- 
dered phase is unstable for all r < c 2 /4a. The quantum- 
critical fluctuations in the disordered phase have been 
discussed by Schmalian and Turlakov. 24 



2. Ordered phase 

In the ordered phase the determination of the Gaussian 
fluctuations is more complicated. Let us parameterize 
the order parameter field as follows, 



M(x) = (m + 6m(x)) 



COs(qz + <j>(x)) 
sin(oz + <p{x)) 
ip(x) 



(3.4a) 



where <fi(x), <p(x), and Sm(x) describe small fluctua- 
tions about the saddle-point solution. Fluctuations of 
the norm of M one expects to be massive, as they are in 
the ferromagnetic case, and an explicit calculation con- 
firms this expectation. We thus can keep the norm of 
M fixed, which means that 5m is quadratic in the small 
fluctuation ip and does not contribute to the Gaussian 
action. In order to treat the <j) and ip fluctuations, it is 
useful to acknowledge that, upon performing a Fourier 
transform, <p(k — 0) corresponds to taking M at k = q, 
while ip and M come at the same wave number. We 
therefore write 



<p(x) = <pi (x) sin qz + Lp 2 (x) cos qz, 



(3.4b) 



where ip\ and Lpi are restricted to containing Fourier com- 
ponents with |fc| <SC q. 25 The Goldstone mode is now ex- 
pected to be a linear combination of <j), ip\ , and ip2 at zero 
wave vector. If we expand Eq. (3.4a) to linear order in 
4> = (fio, substitute this in Eq. (2.1), neglect rapidly fluc- 
tuating Fourier components proportional to e mqz with 
n > 2, and use the equation of motion (3.1a), we obtain 
a Gaussian action 



S (2) N = ^£ E <Pi(pHj(P)<Pj(-P) (3-5a) 

p i=0,l,2 



with a matrix 



7(P) 




-iqpy 
q 2 +p 2 /2 



-iqPx 
q 2 +p 2 /2 



(3.5b) 



The corresponding eigenvalue equation reads 

(H-p 2 )(q 2 +p 2 /2-rf + q 2 pl{q 2 +p 2 /2-ri 

+ q 2 p 2 {p 2 -n) =0. (3.6) 

We see that at p = there is one eigenvalue \i\ = and 
a doubly degenerate eigenvalue ^2,3 — q 2 ■ As expected, 
there thus is one soft (Goldstone) mode in the ordered 
phase. The behavior at nonzero wave vector is easily 
determined by solving Eq. (3.6) perturbatively. The de- 
generacy of fi 2 an d /U3 is lifted, fi 2 ,3 {p — ► 0) = q 2 ± qp z , 
and for [i\ we find 

Mi(P - 0) = p 2 z + pi/2q 2 + 0(p 2 pl). (3.7a) 

The corresponding eigenvector is 

«i(p) = 0(p) - »(P»/?)[l + 0(Pi)]^i(p) 

- l{p x lq)[l + 0{pl)]^ 2 {p).{Z.7h) 
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It has the property 

(«i(p)«i(-p)) = l/am^i(p). (3.7c) 

A comparison with Eq. (2.5) shows that the effective 
soft-mode action has indeed the form that was expected 
from t he an alogy with chiral liquid crystals. If we iden- 
tify y/ a m§ Vi (x) with the generalized phase u(x), the 
coupling constants in Eq. (2.5) are c z = 1 and c± = 1/2. 
Repeating the calculation in the presence of a term that 
breaks the rotational symmetry, e.g., Eq. (2.21), yields a 
result consistent with Eq. (2.22) or (2.23), with b = 0(1). 

IV. NATURE OF THE GOLDSTONE MODE IN 
QUANTUM CHIRAL MAGNETS 

We now turn to the quantum case. Our objective 
is to develop an effective theory for itinerant helimag- 
nets that is analogous to Hertz's theory for itinerant 
ferromagnets. 26 That is, starting from a microscopic 
fcrmionic action we derive a quantum mechanical gener- 
alization of the classical Ginzburg-Landau theory stud- 
ied in the preceding section. The coefficients of this ef- 
fective quantum theory will be given in terms of elec- 
tronic correlation functions, which allows for a semi- 
quantitative analysis of the results. In addition, it pro- 
vides the building blocks for a treatment of quantum he- 
limagnets by means of many-body perturbation theory, 
which will allow to go beyond the treatment at a saddle- 
point/Gaussian level employed in the present paper. Wc 
will show that this theory has a helical ground state given 
by Eqs. (2.2), and consider fluctuations about this state 
to find the Goldstone modes. 



A. Effective action for an itinerant quantum chiral 
magnet 

Consider a partition function 

Z = Jd[4>^} e s ^l (4.1a) 

given by an electronic action of the form 

S$,1>]=S $,il>]+S& t . (4.1b) 

Here Sf nt describes the spin-triplet interaction. S^.xp], 
which we will explicitly specify later, contains all other 
parts of the action, and the action is a functional of 
fcrmionic (i.e., Grassmann-valued) fields ip and ip. The 
spin-triplet interaction we take to have the form 

1 f f 1/T 

S L = 2 / dx d y / dT n K x > T ) A ^ ( x - y) n l (f i T )- 

(4.2a) 

Here and in what follows summation over repeated spin 
indices is implied, x and y denote the position in real 



space, r is the imaginary time variable, and n l s (x,T) de- 
note the components of the electronic spin-density field 
n s (x, t). The interaction amplitude A is given by 

Aij(x -y) = 5ijT t 5(x - y) + e l]k C k (x -y). (4.2b) 

The first term, with a point-like amplitude I\, is the 
usual Hubbard interaction. The second term involves a 
cross-product n s (x) x n s (y) and cannot exist in a ho- 
mogeneous electron system, which in particular is invari- 
ant under spatial inversions. Dzyaloshinski 4 and Moriya 5 
have shown that such a term arises from the spin-orbit in- 
teraction in lattices that lack inversion symmetry. After 
coarse graining, it will then also be present in an effective 
continuum theory valid at length scales large compared 
to the lattice spacing. In such an effective theory the 
vector C(x — y) is conveniently expanded in powers of 
gradients. The lowest-order term in the gradient expan- 
sion is 

C(x-y) =cT t S{x-y) V + 0(V 2 ), (4.2c) 

with c a constant. The ferromagnetic case 26 can be re- 
covered by putting c = 0. We now perform a Hubbard- 
Stratonovich transformation to decouple the spin-triplet 
interaction. To linear order in the gradients the inverse 
of the matrix A has the same form as A itself, viz., 

1 c 
A Tj ( x ~y) = s ij 7T 5{x-y)-e ijk — 5{x-y) <9 fe +0(V 2 ), 
It It 

(4.2d) 

with dk — d/dxk a spatial derivative. The Hubbard- 
Stratonovich transformation thus produces all of the 
terms one gets in the ferromagnetic case, and in addi- 
tion a term 

~ cT t J dx M(x, t) • (V x M(x, t)) + 0(V 2 ), (4.3) 

where M is the Hubbard-Stratonovich field whose expec- 
tation value is proportional to the magnetization. The 
partition function can then be written in the following 
form. 

Z = J D[$,il>] e So ®M J D[M] e -(rt/2)/<faM(x).Af(x) 

xe -c(r t /2)J diM(x)-(VxM(i)) e r t J dx M(x)-n s (x) 

(4.4) 

Here we have adopted a four-vector notation, x = (x,t) 

and J dx = J dx J^ T dr. 

Now we consider the ordered phase and write 

M(x) = M sp {x) + 5M(x), (4.5) 

with M sp given by Eq. (2.2b). The parameters mo and q 
which characterize M sp will still have to be determined. 
By substituting Eq. (4.5) in Eq. (4.4) and formally inte- 
grating out the fermions we can write the partition func- 
tion 




(4.6a) 
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with A an effective action for the order-parameter fluc- 
tuations, 

A[SM] = -InZo + y J dx M(x) ■ M(x) 

+ ~Y J dx M ^ ' *- v x M ( x ^ 



_ ln ( e r t j dx SM( x y na (x)\ (46b) 

/ So 



Here 



SolV*,^] = SotyM +T t J dx M sp (x) ■ n s (x) (4.7a) 

is a reference ensemble action for electrons described by 
So in an effective external magnetic field 



H(x)=T t M sp (x) 



(4.7b) 



Only the Zeeman term due to the effective external field 
is included in the reference ensemble. Zo is the partition 
function of the reference ensemble, 

Z = J e s °®M, (4.7c) 

and (. . .)s denotes an average with respect to the action 
So- 

The effective action A can be expanded in a Landau ex- 
pansion in powers of SM. To quadratic order this yields 



A[5M] = J dx Yf\x)6Mi{x) 

+ ±Jdxdy 5M l (x)T ( i f(x,y)5M 3 (y) 
+ 0(SM 3 ). (4.8a) 

with vertices 

rfOc) =r t (l-cg)M* p (a;)-r t «(*)) So (4.8b) 
and 

r\f(x,y) = 5ij 5(x - y)T t - e ijk S(x - y)T t cd k 

-X^{x,y)Tl (4.8c) 

Here 



Xo(x,y) = (nl{x)nl(y))] 



(4.8d) 



is the spin susceptibility in the reference ensemble. The 
superscript c in Eq. (4.8d) indicates that only connected 
diagrams contribute to this correlation function. 



B. Properties of the reference ensemble 

In order for the formalism developed in the previous 
subsection to be useful we need to determine the proper- 
ties of the reference ensemble. As we will see in Sec. IV E, 



and in a forthcoming paper, 12 the reference ensemble is 
not only necessary for the present formal developments, 
but also forms the basis for calculating all of the thermo- 
dynamic and transport properties of helimagnets. This 
is because the reference ensemble, rather than just being 
a useful artifact, has a precise physical interpretation: it 
incorporates long-range helical order in a fermionic ac- 
tion at a mean-field level. 

We first need to specify the action So. For simplicity, 
we neglect the spin-singlet interaction contained in So 
and consider free electrons with a Green function 



G {k, iu> n ) = 1 1 (ito n - £fc) 



(4.9a) 



Here co n = 2nT(n + 1/2) is a fermionic Matsubara fre- 
quency, and 



£ fc = fc 2 /2m c - e F 



(4.9b) 



with m e the effective mass of the electrons and £f the 
chemical potential or Fermi energy. (Here we neglect 
spin-orbit interaction effects discussed in Ref. 27 as well 
as quenched disorder.) For later reference we also define 
the Fermi wave number fcp = y/2m e eF, the Fermi velocity 
vf = kp/mc, and the density of states per spin on the 
Fermi surface N F = fc F m e /27r 2 in the ensemble S - 



1. Equation of state 

The equation of state can be determined from the 
requirement 28 



(SM(x)) = 0, 



(4.10a) 



where (...) denotes an average with respect to the effec- 
tive action A. To zero-loop order, this condition reads 



(1 - cq) M sp {x) - {n s (x)) So = 0. 



(4.10b) 



The zero-loop order or mean-field equation of state is thus 
determined by the magnetization of the reference ensem- 
ble induced by the effective external field T t M sp (x), Eq. 
(4.7b). The latter is given by the effective field times a 
generalized Lindhardt function. The result is 



1 — cq 



n -i 



G 1 (p, iui n ) G 1 (p - q, iuj n ) - A 2 

(4.10c) 

where A = mo T t is the exchange splitting or Stoner gap. 
Notice that this provides only one relation between A 
and the pitch wave number q. The latter still has to be 
determined from minimizing the free energy, as in the 
classical case, Sec. Ill A. 
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2. Green function 



4- Ferromagnetic limit 



The basic building block for correlation functions of the 
reference ensemble is the Green function associated with 
the action 5 , Eq. (4.7a). With S as specified above, 
the latter reads explicitly 

S a \i>M = j dxdy$(x)G-\x,y)Tl>{y). (4.11) 
with an inverse Green function 



+r t M sp (:E) • er 



5(x-y). (4.12a) 



Here cr = (a\, <T2, as) denotes the Pauli matrices, and cto 
is the 2x2 unit matrix. Upon Fourier transformation we 
have 

G kp( iuj ri) = S k , P GQ 1 {k,ii J j n )(7( S +T t M sp (fc -p) • cr. 

(4.12b) 

The result of the inversion problem is 

G k p{iu n ) = 5kp[<T+-a + (k,q;iu> n ) +a- + a-(k,q;iu> n )] 
+S k+q ,p cr + b+(k, q; iuj n ) 
+Sk-q,p&-b-{k,q;iu> n ), (4.13a) 

where 

Gq x (k ± q,iu) n ) 



a±(k,q;iu n ) = 
b±(k,q;iu>„) 



G 1 (k,iuj n )G 1 (k±q,iuj n ) - A 2 

(4.13b) 

—A 

G " 1 (fc,zw n )Go 1 (fc±q,w„) - A 2 

(4.13c) 

Here <r± = (<7i±i<72)/2, <r^ = er + (7_, and cr |_ = cr_(T + . 

3. Spin susceptibility 



It is illustrative to check the ferromagnetic limit, q — > 
0, at this point. In this case Af sp = (mo, 0,0) becomes 
position independent, and both the Green function and 
the reference ensemble spin susceptibility become diag- 
onal in momentum space. For zero momentum and fre- 
quency, the latter is also diagonal in spin space, 

Xo!q=o( fe 'P;^n) = Sk,pXo, q = ( k ^^n), (4.15a) 

Xo%=o(0, i0) - % [6n xl + (1 - S a ) X t] ■ (4.15b) 

The static and homogeneous transverse susceptibility \t 
of the reference ensemble is related to the magnetization 
by a Ward identity 28 ' 29 (remember that T t M sp is the 
effective field in the reference ensemble, see Eq. (4.7b)) 



(n s (x)) = r t M sp XT- 



(4.16) 



A calculation of xt by evaluating Eq. (4.14b) for q = 
shows that Eq. (4.16) is the equation of state, Eq. (4.10c), 
for q = 0. Equation (4.10c) thus represents the general- 
ization of this Ward identity to the helimagnetic case. 



C. Gaussian fluctuations I: k — q modes 

We are now in a position to explicitly write down the 
fluctuation action given by Eqs. (4.8). From both the 
phenomenological arguments in Sec. II A and the classi- 
cal field theory in Sec. Ill A we expect the static behavior 
to be correctly described by the fluctuations with wave 
numbers close to the pitch wave number q, while Sec. II B 
suggests that treating the dynamics correctly requires to 
also take into account fluctuations with wave numbers 
near zero. For the sake of transparency we first concen- 
trate on the k = q modes. We will later expand our set 
of modes to study the effects of the k = modes on the 
dynamics. 



Since the reference ensemble describes noninteracting 
electrons, the reference ensemble spin susceptibility fac- 
torizes into a product of two Green functions. Applying 
Wick's theorem to Eq. (4.8d) one obtains 



Xo V) = - tr i a i G ( x > V) °i G (V' x )) . 
or, after a Fourier transform, 



(4.14a) 



-1 



Xo(k,p;iCl n ) = — ^2 T ^2 tr {vi G k',p'(iv n ) 

k' .p' ito n 

xctj G p > +Pik > +k (iu>n + i&n))- (4.14b) 

Here the trace is over the spin degrees of freedom, and 
Q n = lixTn is a bosonic Matsubara frequency. 

From the structure of the Green function, Eq. (4.13a), 
it is obvious that xo is nonzero if k and p differ by zero, 
±<7, or ±2q. The full expression in terms of a± and b± 
is lengthy and given in Appendix B. 



1. Gaussian action 

We parameterize the fluctuations of the order param- 
eter as in the classical case, Eqs. (3.4), but now allow for 
the fields 0, ipi, and ip2 to depend on imaginary time or 
Matsubara frequency To linear order in the fluctuations 
we have 

(— 4>{x) sin(q • x) 
4>(x) cos(q • x) 
ipi(x) sin(q • x) + ipz( x ) cos(q • x) 

(4.17) 

As in the classical case we have anticipated that fluctu- 
ations of the norm of the order parameter are massive. 
The term linear in 5M vanishes due to the saddle-point 
condition, and the Gaussian term can be expressed in 
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terms of integrals by using Eqs. (4.14b) and (B1,B2). 
Using the notation <p(x) = <po(x) as in the classical case, 
we find a Gaussian action 30 

Ai2) i^ = yEE E Vi(P,Mn)ll?( P ,inn) 
p iQ. n i=0,l,2 

X(pj(—p,—iCl n ). (4.18a) 
I 



Here the matrix is the quantum mechanical ana- 
log of Eq. (3.5b), which couples the phase or k = q 
modes among each other. In a four-vector notation, 
k = (k, iSl n ) = (k x ,k y , k z ,iQ, n ) it is given by 



Here 



(1 - cq)/T t - f H {k) -ick y /2T t -ick x /2T t 
7 (,) (*)=| ick y /2T t l/2r t -/ n (fc) -/i 2 (fe) | (4.18b) 

ick x /2T t f 12 (k) l/2r t -/ n (fc) 



Utity = <P<t,<t>{k) +<P</>4>(-k), (4.18c) 
fu(k) = <pu(k) + <pu(-k), (4.18d) 
f 12 (k) = ifa^k) - <p u (-k)], (4.18e) 



where 



G 1 (P- fc i iu m - itt n )G 1 (p-q, iu m ) - A 2 



iui m - iO„)G 1 (p-k- q,iu m - ifi n ) - A 2 ] [G ^p,^) G *(p- q, iui m ) - A 2 ] ' 

(4.18f) 

^ n ( fc ) = — ( G o 1 ip~ k -, i ^m - ifl, l )GQ 1 (p + q,iuj m ) + gp^jg - fc q,^m - ^n)Gp ^P,^) - 2A 2 

4 7 p [Go 1 (p-fc,ia; m -ifi„)G'o 1 (p-fc-q,w m -2f7„)-A 2 ] [G^ *(p, zw m ) G^p + q, iw m ) - A 2 ] ' 

(4.18g) 



withj^(lA0£ p T£^. 

In contrast to the classical case, here it is not obvious 
that the Gaussian vertex, Eq. (4.18b) has a zero eigen- 
value. To see that it does, we invoke the equation of state 
(4.10c). By comparing this with Eqs. (4.18c, 4.18f), we 
see that 



C, yields 

*n„) - /,^(o,io)-2Jv F 



a 



AY 

2k F J 



-R ( 

'' Ue F 



1<P 



l"n| kl 

4e F 2fc F |fc| 



(4.20a) 



Similarly, 



l-c«-r t /^(0,i0) = 0. (4.19a) 



l/2-r t /n(0,»0) = cg/2. (4.19b) 



Since c oc q (see Eqs. (3.2a) and (4.21) below), it fol- 
lows that the quantum mechanical vertex "f^ q \p, id n ) 
has the same structure as its classical counterpart, Eq. 
(3.5b), except for an additional frequency dependence in 
the quantum mechanical case. 

To determine the eigenvalues we need to evaluate the 
integrals to lowest nontrivial order in the wave vector and 
the frequency. A complete calculation is rather difficult, 
and we restrict ourselves to the limit A 3> qvp. The 
calculation, the details of which we relegate to Appendix 



fu{k,m n ) = /n(0,tO)-JV F 

2 



2fa 



\4e F 



■7i 



k\ 



fi2(k,iQ n ) 
Here 



-i2N F a 



4e F (2fc F ) 2 

qk z 



(2fc F ) 2 



(4.20b) 
(4.20c) 



a = 1/3, (4.20d) 

(3 = 4e|/A 2 , (4.20e) 

70 = 7T(qv F ) 2 /8X 2 , (4.20f) 

71 = -4 70 (fc F / g ) 3 . (4.20g) 

These expressions are valid for \Q n \ <C A <C e F , \k\ <C 

ij « tp, and qvp <C A. The damping terms have the 
form shown if, in addition, IOI < j^lfel. 
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We now also can express the pitch wave number q in 
terms of the parameters of our model. The minimization 
of the saddle-point free energy proceeds as in the classical 
case, and by comparing Eq. (3.2a) with Eqs. (4.17,4.18a, 
4.18b, 4.19a, 4.20a), we find 

q = ck F /N F T t a. (4.21) 

2. Eigenvalue problem 

The diagonalization of the matrix Y q \ Eq. (4.18b), 
can be done perturbatively as in the classical case. The 



soft (Goldstone) mode, i.e., the eigenvector correspond- 
ing to the smallest eigenvalue, is 

v(k,iCl n ) = c/)(k,in n ) - i(k y /q)[l + 0(k 2 ± )](f 1 (k,in n ) 
-i(k x /q)[l + O(fci)] P2(fc,in„). (4.22) 

The helimagnon is proportional to the u-w-correlation 
function, which in turn is proportional to the inverse of 
the smallest eigenvalue of the matrix 7^. With k = 
^/ak/2k F , Q = ^/aq/2k F , lu = ^J]3Vt n /Ae F = fi„/2A, 
c = l4>/V a Pi and c i — li/ a VP, the latter reads 



7 (q) (fc) = 2N F 



' + C0|CJ| 

i Q Ky 



K 



Q 2 + 



K 2 + 



i Q Ky 

-iQn z 



Q 2 + 



■k 2 + 



,2 



\ci\u\ 



(4.23) 



The eigenvalue equation, which is the quantum me- 
chanical generalization of Eq. (3.6), can be simplified if 
we anticipate, from Sees. II and IIIB, that the smallest 
eigenvalue scales as fi ~ k 2 ~ k\ <~ uj 2 . Keeping only 
terms up to 0(k 2 ), the eigenvalue equation reads 

^-k 2 -u> 2 - C4> nj£) (q 4 + qV ± ) 

+ Q 2 nl(Q 2 + ?pj =0. (4.24) 

This result has several interesting aspects, which will 
be very useful when we generalize the theory to include 
the k = modes in Sec. IV D. First, of the function / n , 
Eq. (4.20b), only the constant and the term proportional 
to k\ ~ k z contribute to the leading terms in the eigen- 
value. In particular, the damping term proportional to 
71 , which has a potential to lead to an overdamping of the 
Goldstone mode, does not contribute. Second, the func- 
tion /12, which describes the coupling between the phase 
modes ipi and (f2, does not contribute to the leading re- 
sult. Indeed, the only role of tpi and </?2 is to subtract the 
k\ contribution from the eigenvalue, and this is accom- 
plished entirely by the purely static matrix elements in 
Eq. (4.18b) which do not depend on reference ensemble 
correlation functions. 

For the smallest eigenvalue we find from Eq. (4.24) 

k 4 k 2 
H{n -» 0, oj -» 0) = k 2 z + uj 2 + ^ + c |w| j-i + 0(4). 

(4-25) 

Since the Goldstone correlation function is propor- 
tional to the inverse of the smallest eigenvalue, this re- 
sult has indeed the functional form we expect from the 
phenomcnological treatment in Sec. II, see Eqs. (2.18a, 
2.19a). At zero frequency it also is consistent with the 



result of the classical field theory in Sec. III. However, in 
contrast to Eq. (2.13) the mass of the zero- wave number 
mode does not enter in Eq. (4.25), which suggests that 
the frequency scale is not correctly described yet. This 
was to be expected, see the remarks at the start of Sec. 
IV C. We will see in Sec. IV D that the k = mode 
has to be included in the analysis to obtain the correct 
prefactor of the uj 2 term in the smallest eigenvalue, in 
agreement with the phenomenological analysis. 

3. Ferromagnetic limit 

Before we expand our set of modes it is again illustra- 
tive to consider the ferromagnetic limit q — 0. Going 
back to Eqs. (4.17), we see that in this limit ipi disap- 
pears, and qb and ip2 play the roles of the two independent 
field components 7Ti and n 2 in a ferromagnetic nonlinear 
sigma model, which are both soft. 31 The Gaussian action 
is given by the reference ensemble spin susceptibility at 
q = 0, and one finds 32 

A {2) [c^, V2 ] = \N F T 2 

p,ifi„ i=l,2 

xttj (4.26) 
with 7Ti = (j) an d T2 = V?2, and a matrix 

%(p,^)=(X) V)- ( 4 - 2? ) 

Comparing with the eigenvalue problem given by Eq. 
(4.23) we see that the latter does not correctly repro- 
duce the ferromagnetic result upon dropping ipi and let- 
ting q — > 0. This is not surprising, since the gradient 
expansion implicit in the effective field theory approach 
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implies that we are restricted to wave numbers small com- 
pared to q. That is, the validity of Eq. (4.23) shrinks to 
zero as q — > 0. However, it raises the following ques- 
tion. The leading frequency dependence of the ferromag- 
netic magnon is determined by the off-diagonal elements 
of the matrix 7. Similarly, the time-dependent Ginzburg- 
Landau theory of Sec. II B suggests that the leading fre- 
quency dependence of the helimagnon is produced by the 
coupling of the phase mode to the homogeneous magne- 
tization. Although the latter is massive at wave number 
q, its conserved character makes it important for the dy- 
namics. In contrast, in the treatment above the leading 
frequency dependence was produced by the phase corre- 
lation functions. We therefore extend our set of modes 
to allow for a coupling between the phase modes, which 
represent k = q fluctuations of the magnetization, and 
the homogeneous magnetization. 



D. Gaussian fluctuations II: coupling between 
k = q modes and k — modes 

The discussion in the preceding subsection suggests to 
generalize the expression for the magnetization fluctua- 
tions, Eq. (4.17), by writing 



SM(x) = m 



I —<f>(x) sin(q • x) 

4>{x) COs(q • x) + TT2(x) 

fi(x) sin(q • x) + ^{x) cos(q • x) 

V +7Tl(x) 



That is, we add the k 
fluctuations. 33 



(4.28) 

fluctuations to the k = q 



1. Gaussian action 

If one repeats the development of the previous subsec- 
tion for the current set of five modes, one finds that of the 
two k = modes only 7Ti contributes to the leading terms 
in the eigenvalue problem. 7r 2 does not couple to <f>, and 



its couplings to ipi and <f2 produce only higher order cor- 
rections. 7Ti, on the other hand, couples to <fi in a way that 
preserves the off-diagonal frequency terms characteristic 
for the ferromagnetic problem, see Eq. (4.27), and needs 
to be kept. These observations are in agreement with 
the phenomenological theory of Sec. II B, where only the 
3-component of the homogeneous magnetization coupled 
to the phase mode. They also are consistent with the 
observation in Sec. IV C 2 that ipi and (f2 serve only to 
provide the correct static structure of the theory. We 
thus drop 7T2 and consider the 4x4 problem given by the 
phase modes plus w±. 

The Gaussian action that generalizes Eqs. (4.18) now 
reads 



p i=0 



xipj(-p,-iQ n ). (4.29a) 



Here (ps = 7Ti, and the matrix -y(9>°) ) which couples the 
phase modes to m, reads, in a block matrix notation, 



7 («.0)(jfe) 



V i/tyi(fc) 



-ih<f,i(k) \ 





i/r t - 5 n(fc) J 



(4.29b) 



with 7^ from Eq. (4.18b). In addition to the functions 
defined in Eqs. (4.18), we need 

gu{k,iSl n ) = 4(^ii (k - q,iSl n ), (4.29c) 

with Lpn from Eq. (4.18g), and 

/i0i(fe,iO„) = rfai(fc,if2„) - ffai(-fc, -i£l n ), (4.29d) 

where 



^i(fc) = A 



G 1 (p- k,iui m - ifi n ) ~G 1 (p- q,iuj m ) 



k, iuj m - iSl n ) G 1 (p- k - q, iw m - iSl n ) - A 2 ] [G 1 (p, iuj m ) G 1 (p- q, iu m ) - A 2 ] 

(4.29e) 

I 



Performing these integrals in the limit A 3> qvp yields the following leading contributions (see Appendix C), 



gn(k,in n ) = l/r t -2AT F 



. (4.30a) 



io. n y 1 \n n \ 2k F k\ 



4e F J 2 19 4e F |fc| 3 



^l(fc,in„) = -2N ¥ yff}—±. 

4ep 



(4.30b) 
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2. Eigenvalue problem 



where 



Now consider the generalization of Eq. (4.23). It is 
obvious from Eqs. (4.29b, 4.30) that the mode m will 
contribute to the eigenvalue a term proportional to SI 2 
whose prefactor is large compared to the one in Eq. (4.25) 
by a factor of (kp/q) 2 . This is because the mass of n\ is 
proportional to (q/k-p) 2 . Consequently, one can neglect 
the terms proportional to u 2 in the diagonal elements of 
7^), Eq. (4.23). The leading damping term, however, 
is still given by Eq. (4.20a). Dropping all terms 
that yield contributions of higher order than k 2 in the 
eigenvalue equation, the matrix Eq. (4.29b), reads, 

in the notation of Sec. IV C 2, 



,(-3.0) 



(k) = 

(n 2 



2N r 



+ c 4>\ UJ \\i\ -iQk v 

iQK>y Q 

—i(iu>) 



2 



iQn z 

f + \ «i 





i(iuj) \ 



Q 2 

(4.31) 



Mk) = A -f- \k 2 /(2k F ) 2 + - k\/(2qk F ) 2 , (4.33b) 



and 



3k F 



1[ ' 3 £F (2fc F ) 4 2fc F |fe| 



(4.33c) 



The smallest eigenvalue is 



This has the same functional form as the result of the 
phenomcnological treatment in Sec. II, sec Eqs. (2.18a, 
2.19a). The current microscopic derivation reveals in ad- 
dition that the prefactor of the damping coefficient is 
smaller than that of the resonance frequency by at least 
a factor of (qv F /X)(q/k F ). The Goldstone mode is thus 
propagating and weakly damped for all orientations of 
the wave vector. 34 As in the classical case, adding cubic 
anisotropic terms, e.g., the quantum mechanical general- 
ization of Eq. (2.21), leads to a soft-mode energy of the 
form of Eq. (2.23). 



The Goldstone mode in the presence of quenched 
disorder 



fi(n — > 0, u — > 0) 



K 2 + 



+ 



+ C0 \uj\ 



Q 2 2Q 2 

and the corresponding eigenvector reads 



K 



-0(4), 

(4.32a) 



v(k,iil n ) = </>(k,ifl n ) - i(Ky/Q)ipi(k,iQ n ) 

-i(n x /Q) <f2{k, iil n ) + (4^/Q 2 ) 7ri(fe, iQ, n ). 

(4.32b) 

Notice that the origin of the leading frequency de- 
pendence is consistent with the ferromagnetic limit, Eq. 
(4.27), and, more importantly, with the time-dependent 
Ginzburg-Landau theory in Sec. II B. Namely, it is pro- 
duced by the coupling of the 3-component of the k = 
magnetization fluctuation to the phase mode at k = q. 
The frequency dependence in Eq. (4.23), on the other 
hand, corresponded to second order time derivative cor- 
rections to the kinetic equation (2.6). 



E. The Goldstone mode, and the spin susceptibility 

1. The Goldstone mode in the clean limit 

We are now in a position to determine the Gold- 
stone mode. Defining the latter as g(k, iQ n ) = 
(V2A^FV / 3fcF / q) v(k, ifl n ), and returning to ordinary 
units, we find from Eqs. (4.32) 

i 

(g(k,ifl n ) g(-k,-iQ. n )) 



-( ^ ^l n ) 2 + LJ 2 (k) + \n n \-f(k) , 



Equation (4.33c) holds for clean systems. As we have 
mentioned in the context of Eqs. (2.19), the structure 
of the damping term is expected to change qualitatively 
in the presence of even weak impurity scattering. Let 
Timp be the inelastic scattering rate due to quenched im- 
purities. Then the bare Green function Go, Eq. (4.9a), 
acquires a finite lifetime, 35 

G (k, iuj n ) = l/(iw n - 6e + (i/2T im p)sgn uj n ). (4.34) 

For weak disorder, eFTi mp 3> 1, the leading effect is that 
the hydrodynamic singularity in the generalized Lind- 
hard function (f>\!, Eq. (Clf), is now protected. The net 



effect is the replacement l/i;F|fc| 



in the damping 



term in /<^, Eq. (4.20a), and hence in the damping co- 
efficient 7(fe). The Goldstone mode is then given by Eq. 
(4.33a), with the resonance frequency still given by Eq. 
(4.33b), and 



li k ) = ~T e F (fFTimp) 



(2fc F ) 4 {2k F f 



(4.35) 



(4.33a) 



This is again consistent with the result of the phenomeno- 
logical treatment in Sec. II, see Eq. (2.19b). Since 
7(fe) ~ k 2 <C Wo(fc) ~ k z , the mode is again propagating. 



3. The physical spin susceptibility 

The helimagnon correlation function is simply related 
to the physical susceptibility, which is directly measur- 
able by inelastic neutron scattering. 36-38 To express the 
latter in terms of order-parameter correlation functions 
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we generalize the partition function, Eq. (4.1a), to a gen- 
crating functional for spin density correlation functions, 



one finds 



Z\j) = J D$,i>] e s ^+! 



dxj(x)-n s (x) 



(4.36) 



where j(x) is a source field. The spin susceptibility is 
given by 



(nl(x)nl(y)) s 
6 2 



(ni(x)) s (ni(y))s 
\nZ[j] 



3=0 



= Tl({A-HM).{x) (A- 1 SM) ] (y)) A 

(4.37) 

Here A' 1 is the matrix given in Eq. (4. 2d), and the last 
average is to be taken with respect to the effective action 
A, Eqs. (4.8). 

The chiral nature of the helix is also reflected in the 
spin fluctuations and will become manifest if the spin 
susceptibility is measured with circularly polarized neu- 
trons. To see this, it is useful to define magnetization 
fluctuations 6M± — 6 Mi ± i<5M 2 . Similarly, we define 
gradient operators d± = d x ± id y . With these definitions 



r t (A- 1 SM) ± (x) 

r t (a- 1 sm) z (x) 



6M ± {x)±icd z 6M ± (x) 
Ticd±5M z (x), (4.38a) 
5M z (x)+i^(d+5M-{x) 

-dSM+(x)). (4.38b) 



By means of Eq. (4.28) we can express the components 
of SM in terms of tpo = </>, <fii, <fi2, and ip% = m, 



5M±(x) 
SM z (x) 



±im (t>(x) e ±iq - x , 
m 7Ti(a;) + [((f2(x 



(4.39a) 
iVi(x))e ±l i- x 
+c.c] , (4.39b) 



where c.c. denotes the complex conjugate of the preced- 
ing expression. The elements of the susceptibility tensor 
can therefore be expressed in terms of the correlation 
functions of the <pi, which we denote by x<f>4> etc The 
various correlation functions can be determined from the 
Gaussian action A^ 2 \ Eq. (4.29a), with the vertex func- 
tion ^(9,0) gj ven by (4.31). The inverse of the latter 
matrix reads 



/ 



( 7 (^°)(p,^„) 



2JV F Q 3 /i(p,ifi„) 



Q iQ Ky 

'%Q Ky Qfcy 



5 n - :C n y 
-iUJKy 



iQ^ k x 

QkI 
—iuiK x Qk 2 



Qui 

— iLUKy 

k, 4 J2Q + Qc^\k 2 z /\k\ J 



(4.40) 



Here we use the same notation as is Eqs. (4.23) and 
(4.31), and corrections to each matrix elements are one 
power higher in k z <~ «j_ ~ u> than the terms shown. All 
of the correlation functions that determine the spin sus- 
ceptibility are proportional to the inverse of the smallest 
eigenvalue fi, Eq. (4.32a), which scales as u 2 . (cfxfi) is the 
softest; it scales as 1/oj 2 . The autocorrelation functions 
of 4>\ and (j>2 have an additional factor of n\ ~ uj in the 
numerator, and thus scale as 1/w, and so does {(j>\<j>2)- 
The autocorrelation function of 7ri scales as a constant. 
The mixed correlations {<jxj>\) and (tjxfe) scale as 1/oj 3 / 2 , 
and the mixed correlations {<fmi} and ((fii^i) scale as 
1/lu and I/uj 1 / 2 , respectively. 

Defining the momentum and frequency dependent spin 
susceptibility by 

r r 1/T 

xi j {k,p;iQ n ) = J dx dy e lk x - lp v J dr e~ innT 

x X i j (x,y-,T), (4.41) 

we find that Xs an d Xs + reflect the strongest hydro- 
dynamic contribution, which is given in terms of the <p-(p 



correlation function at wave vector ±q, 

xt T (k,p;itt n ) = Sk P mlx4>4>( k ± Q,i^n) 

+ (less leading terms). (4.42a) 

The "less leading terms" in Eq. (4.42a) reflect the terms 
proportional to c in Eq. (4.38a). They are either less 
singular than the leading term, which scales as 1/oj 2 , or 
have a prefactor that is small by a factor of cq oc q 2 . xt + 
and Xs also are proportional to x<f>4>, but they arc not 
diagonal in the momenta, 

X^Cfc.PS^n) = -5p,k±2qmlx<t>4>{k ± q,iQ. n ). (4.42b) 

Only the terms with k — p contribute to the neutron 
scattering cross section. Right and left circularly polar- 
ized neutrons will therefore see the hydrodynamic sin- 
gularity only at wave vector k = q and k = —q, respec- 
tively. Unpolarized neutrons will see symmetric contribu- 
tions at k = ±q. For instance, the xx component of the 
susceptibility tensor, xl 1 = (xt ~ + X7 + + xt + + X7~)/ 4 
is given by 
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(4.43) 



with the first term in the square brackets contributing to 
the neutron scattering cross section. 

The longitudinal component x 33 is l ess singular than 
Xs since the contribution of <j> to (A~ 1 5M) Z is sup- 
pressed by a transverse gradient, see Eqs. (4.38b, 4.39a). 
Accordingly, the leading hydrodynamic contribution to 
Xs 3 at k = ±q scales as (<pi<pi) ~ 1/uj. There is also a 
weak signature of the Goldstone mode in the vicinity of 
k = due to the contribution of (ttiiti) to x 33 , but this 
scales only as ui°. 

All of these results are in agreement with what one 
expects from Sec. II D. 

V. DISCUSSION AND CONCLUSION 

Let us summarize our results. We have developed a 
framework for a theoretical description of itinerant quan- 
tum helimagnets that is analogous to Hertz's treatment 
of ferromagnets. 26 We have analyzed this theory in the 
helically ordered phase at a mean-field/Gaussian level 
analogous to Stoner theory. As in the ferromagnetic case 
there is a split Fermi surface, with the splitting propor- 
tional to the amplitude of the helically modulated mag- 
netization. We then focused on the Goldstone mode, 
or helimagnon, that results from the spontaneous break- 
ing of the translational symmetry in the helical phase. 
We have found that the helimagnon is a propagating, 
weakly damped mode with a strongly anisotropic dis- 
persion relation. The frequency scales linearly with the 
wave number for wave vectors parallel to the pitch of the 
helix, but quadratically for wave vectors in the trans- 
verse direction. In this sense the helimagnon behaves 
like an antifcrromagnctic magnon in the longitudinal di- 
rection, but like a ferromagnetic one in the transverse di- 
rection. This anisotropy is analogous to the situation in 
chiral liquid crystals, and indeed the results of our micro- 
scopic theory are qualitatively reproduced by combining 
an educated guess of the statics, inferred from the liquid- 
crystal case, with a phenomenological time-dependent 
Ginzburg-Landau theory for the dynamics. The struc- 
ture of the microscopic field theory is in one-to-one cor- 
respondence with the structure of the phenomenological 
theory. Particle-hole excitations provide a damping of 
the helimagnon, with a damping coefficient proportional 
to the pitch wave number to the fourth power. 

Our continuum theory ignores the spin-orbit coupling 
of the electron spins to the underlying lattice which, in 
conjunction with crystal-field effects, will change the dis- 
persion relation of the Goldstone mode at very small 
wave numbers or frequencies. In contrast to the case 
of ferromagnetic and antifcrromagnetic magnons, how- 



ever, breaking the spin rotational symmetry does not give 
the helimagnons a mass; it just changes the dispersion 
relation at unobservably small wave numbers, making 
it less soft. Consistent with this, soft helimagnons are 
observable by neutron scattering, 38 as are ferromagnetic 
and antifcrromagnetic magnons. The neutron scattering 
cross section is proportional to the magnetic structure 
factor, which in turn is simply related to the spin suscep- 
tibility, which we have shown to be proportional to the 
helimagnon correlation function. 

Let us now discuss some observable properties, us- 
ing parameter values appropriate for MnSi. MnSi has a 
Fermi temperature Tp w 23, 200 K, 39 a pitch wave num- 
ber q w 0.035 A , 40 and an effective electron mass, av- 
eraged over the Fermi surface, m rts 4m , 41 with to 
the free electron mass. In a nearly-free-electron model 

this leads to k F w 1.45 A - and qv F /k B « 1,000K. The 
value of the exchange splitting A is less clear. The large 
ordered moment of about 0.4 /zb per formula unit 11 sug- 
gests an exchange splitting that is a substantial fraction 
of Tp. This is hard to reconcile with the low Curie tem- 
perature Tq — 29.5 K at ambient pressure. Ref. 41 found 
an exchange splitting A/fee ~ 520 K, which is hard to rec- 
oncile with the large ordered moment. We recall that the 
helical order is caused by the (weak) spin-orbit interac- 
tion, which suggests a clear separation of energy scales, 
and in particular A 3> qvp. In judging such estimates 
one needs to keep in mind that MnSi is a fairly strong 
magnet (as evidenced by the large ordered moment) with 
a complicated Fermi surface with both electron and hole 
orbits, and a nearly-free-clectron model as well as a weak- 
coupling Stoner picture are of limited applicability. 

One question that arises is the value of qvp/X ap- 
propriate for MnSi (recall that our calculation is for 
qvp/X -C 1). This hinges on the value of A. If one ac- 
cepts a sizeable value of A (in units of the Fermi tempera- 
ture) as suggested by the large value of the ordered mag- 
netic moment, then qv-p/X <C 1. If one accepts the much 
smaller value for A obtained in Ref. 41, then qvp/\ w 1. 
Even in the latter case we expect our considerations to 
still apply qualitatively, although the lack of a clear sepa- 
ration of energy scales would make a quantitative analysis 
more difficult. 

The anisotropic helimagnon is well defined for wave 
numbers \k\ < q. In MnSi, q w 0.035 A , which is 
accessible by inelastic neutron scattering. 11 An estimate 
of the helimagnon excitation energy in MnSi at this wave 
number can be obtained from Ref. 42. In this experiment, 
a small magnetic field was used to destroy the helix. The 
resulting ferromagnetic magnons at a wave number equal 
to q had an energy of about 300 mK above the field- 
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induced gap. An estimate of the helimagnon energy at 
the same wave number from Eq. (4.33b) yields the same 
order of magnitude if A is a substantial fraction of the 
Fermi energy. The prediction is thus for an anisotropic 
helimagnon, which at a wave number on the order of q, 
and an energy on the order of 300 mK, will cross over to 
a ferromagnetic magnon. The damping is expected to be 
weak, especially in systems with some quenched disorder, 
although in ultraclean systems it may be strongly wave 
vector dependent. 34 

One expects the low-energy helimagnon mode to have 
an appreciable effect on other observables such as the 
specific heat, or the electrical resistivity. This is indeed 
the case, and we will discuss these effects in a separate 
paper. 12 

In conclusion, the present paper provides a general 
many-body formalism for itinerant quantum helimag- 
nets, which can be used for calculating any observable 
of interest, with the single-particle Green function and 
the helimagnon propagator as building blocks. Observ- 
ables of obvious interest include the specific heat, the 
quasi-particle relaxation time, and the resistivity. Cal- 
culations of these quantities within the framework of the 
present theory will be reported in a separate paper. 12 
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APPENDIX A: 
ROTATIONAL-INVARIANCE-BREAKING 
TERMS IN THE GINZBURG LANDAU 
EXPANSION 



In addition to the term considered in Eq. (2.21), there 
are other terms that break the rotational symmetry and 
contribute to the same order in the spin-orbit interaction 
strength 5 so- The detailed structure of such terms in the 
action depends on the precise lattice structure. For defi- 
niteness we assume the cubic P2i3 space group realized, 
for instance, in MnSi. 

The leading terms which break the rotational symme- 
try are of the form 6,7 



SS = j dx {oi [(d 2 M(x)) 2 + (d 2 M{x)) 2 + (d 2 M(x)) 2 ] + g 2 so [a 2 (d x M y (x)f + a 3 (d x M z (x)) 2 + cycl.] 

[(M x (x)f + {M v (x)f + (M z {x)f}). (Al) 



3so a 4 



Here ai, a 2 , a 3 and 014 are constants which remain finite 
for vanishing spin-orbit coupling, <?so ~~ * 0, and cycl. 
denotes cyclic permutations of x, y, and z. Omitted 
terms like (MjM 2 + cycl.) or \{d x M x ) 2 + cycl.] can be 
obtained by adding rotationally invariant terms to SS. 
We note that the Dzyaloshinski-Moriya interaction, i.e. 
the constant c in Eq. (2.1), the pitch wave number q, and 
therefore all typical momenta, are linear in g so ■ Conse- 
quently, all terms in Eq. (Al) contribute only to order 
g g . They are therefore small compared to the energy 
gain obtained by forming the helical state, which is is of 
order q 2 oc g% Q . 

The main effect of SS is to pin the direction of the 
spiral to some high-symmetry direction, namely, either 



(1,1,1) and equivalent directions, or (1,0,0) and equiva- 
lent directions. 6,7 In addition, they change the dispersion 
relation of the helical Goldstone mode at extremely small 
wave numbers. For one particular term this has been 
demonstrated in Sec. HE; all other terms have qualita- 
tively the same effect. Notice that in a /erromagnet the 
effect on the Goldstone mode is stronger. For instance, 
the cubic anisotropy (the term with coupling constant 
a 4 in Eq. (Al)) gives the ferromagnetic magnons a true 
mass, leaving no soft modes. This is a result of the fact 
that in a ferromagnetic state the translational invariancc 
is not spontaneously broken, while in a helimagnetic one 
it is. 
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APPENDIX B: THE REFERENCE ENSEMBLE SPIN SUSCEPTIBILITY 

Substituting Eq. (4.13a) in Eq. (4.14b) and performing the spin traces yields 

xl 3 (k,p;iSl n ) = -rrE T Er fc ^ ( S <^ + )*j ■ a +( k ' ~ fc > «?; - i^n) a+(k', q; iu n ) 



+S k - 



+ 

+(^ + 

■P (^ab + 



.. a+(k' - fe,q;iw : 
i:j a-{k' -k,q 
.. a-(k' -k,q 
ij b+(k' -k,q 
..b.{k'-k,q 
ij a+(k' k,q 

.,b+(k' -k,q 

6 J 

ij M fe ' k ,Q 
..a+{k'-k,q 

ij a -( k ' k ,q 

..b-{k'-k,q 



b_(k'-k,q 



lUJr, 

iuj n 
iu)„ 
iuj„ 
iu> r , 

iuj„ 
iu>„ 
iu> r , 
iuj n 
iuj„ 

iu>„ 



■ iCl n )a-(k' ,q;iw n ) 

- in n )a + (k' ,q;iw n ) 

- ifl n )a-(k',q;iw n ) 
ifl n )b-(k' + q,q;iuj n ) 

■iCl n )b+(k' - q,q;iw n ) 

- iQ n ) b + (k' - q, q; iuj n ) 

- iCl n )b+(k' - q,q;ico n ) 
iQ n ) a+(fe', q; iuj n ) 

iCl n )a-(k' ,q;iw n ) 

- iCl n )b-(k' + q,q;iu> n ) 

- ifl n )b-(k' + q,q;iu> n ) 
iil n )a + (k / ,q;iuj n ) 

- iVL n )a_(k' ,q;iu) n ) 



+$k-2q, P (^bb + )ij M fc ' - k,q;iu m - ifl n ) b+(k' - q,q;iw n ) 
+5 k +2q., P (^bb~)ij b-{k' - k,q;iw m - ifl n )b-(k' + q,q;iu n ) j. 



(Bl) 



The £ symbols denote traces of Pauli matrices 





(£++) =tv(a i a + a_a j a + a_)= I 

1 1 



(B2a) 



{^aa )ij =tl(ai(T + <T-<Tj<T-a + ) = [ i 1 ) , 



(B2b) 



1 i 

( S aa + ) y = tr (ffj (T_ 0+ Oj- + (T_) = I -* 1 I , 









(£"-) =tv(a i a-a + a j a-a + ) = I 



1 



(B2c) 



(B2d) 



(Kb )ij = tr ((Ti a+ CTj ct_) = 








(S» + ) i ,=tr(^a_ - J - -+)= 

1 -1 



1 

( S £% =tr(<j l a + a-a J a+) = I i 

1 



(B2e) 



(B2f) 



(B2g) 








(Ka + ) l3 = tr (a i( r + a ja + a_) = 0, 

1 1 i ' 



(B2h) 
(B2i) 
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APPENDIX C: THE FUNCTIONS fa, fa, fa, gu 

=tv(a i a + a j a-a+) = ( -i | , (B2j) AND ft^ 



In this appendix we show how to evaluate the functions 

f<t>4>, hi, /12, and h^i defined in Eqs. (4.18c - 4.18g) 

/ n / " " " \ /T , 01 s and (4.30) in the limit of lone wavelengths and small 

(Kb ) ii =tr(a i a + a.a j a.)= \ |, (B2k) 




(S- 6 -)..=tr(a iff _ - + ^a_)= [ /], (B21) The function 

/ 

We start with the expression for ip^ in terms of Green 
functions given in Eq. (4.18f). By symmetrizing the de- 
1 \ pendence on k and q, performing the sum over Matsub- 

(^bo + )i, ; = tr (°» a - a j a + a -) = ( — i J , (B2m) ara frequencies, and doing a partial fraction decomposi- 

J tion, ip^ can be written 

o o o\ ^0 (fc,^„) = tp { +J(k,in n ) + tp { ^{-k,-m n ) 

(Sr - ) . . = tr (ai a- <r, ; <r_ a+) = I 0, (B2n) , (-)/, . n n (-)/ , n 

v ba J v v 3 + > I _i j o / ^ ^ k ' iQn > ~ v <t"t> (- fc >-^«). 

(Cla) 



1 i 0\ From Eq. (4.18c) we have 

( s » + )« =tr( f 7 i a + a J <7 + ) = I i -1 , (B2o) 







1 -i 



fa{k,iQ„) = 2 ^(fc.ifin) + <^ ) (-fc,-iQ„) . 

(Clb) 

That is, is given by the symmetric part of (ys^^ alone. 



)a =tr(cj i a-cj j a-) = -i -1 , (B2p) The antisymmetric nart 







The antisymmetric part does not contribute, but we 



list it here for completeness: 



2m c w + W- 



+ 



{q-P-)w+ + {q-p+)w- 



iQ n — k ■ p/m c + w + — w_ iQ n — k ■ p/m c — w + + w 



iCt n — k ■ p/m c — w + — w- iCt n — k ■ p/m c + w + + w- 



2m c w + w- 

The symmetric part, </A+), consists of two parts that are structurally distinct, 



where 



(i) 



(k,ifl n ) = —Y, 



1 \ ^ + A 2 + g(p) 



W + W- 



+ - 



/(£+«>-) 



iQ n — k ■ p/m c — w + — w- iQ n — k ■ p/m c + w+ + w 



. (Clc) 



(Cld) 



(Cle) 



(2) 



+ ■ 



ifl n — k ■ p/m c + w + — w_ ifi„ — k ■ p/m c — w + + W- 
I 



In these expressions £± = £ p± + q 2 /8m e , p± = p± fe/2, and 
and w± = w(p±) with 



g(p) 



1 I. x2 1 



(Clg) 



(Clf) 



(Clh) 



19 



y>00 and ip^ are both generalized Lindhard functions. 
For A = they combine to form a Lindhard function at 
wave vector k + q, while for q = ip^ vanishes and ip^ 
turns into the function that determines the ferromagnetic 
magnon 32 . A structural difference between them is that 
in the hydrodynamic singularity at k = ifl n = that 
is characteristic for the Lindhard function is protected by 
A, while in this is not the case. 

The remaining wave vector integral is difficult, and we 
evaluate it only in the limit A ^> qvp. The Fermi func- 
tions in Eqs. (Cle,Clf) pin the integration momentum p 
to a shifted Fermi surface. For A > qvp one can therefore 
perform a straightforward expansion of the integrand in 
powers of qvp/X. For ipW the leading term is the q = 
contribution. Specifically, in this limit 



(w + w- + A 2 + g(p))/w + W- 



(C2) 



and the calculation reduces to the ferromagnetic case. 
For k <C fc F and \Q n \ <C A, the result is 



¥$(0,iO)+JV F 



k 

2k F 



+ 



~2X 



(C3) 



For (p^ the leading term is of 0(q 2 ), since 

(w + w- - A 2 - g{p))/w+w- -» (g ■ fe) 2 /8m 2 A 2 , (C4) 

and the integral is reduces to a Lindhard function. 
The non-hydrodynamic part provides only corrections of 
0((qvF/X) 2 ) to (p^ , but the hydrodynamic part is qual- 
itatively new and must be kept. We find 

which is valid for \Q n \ -C vpk. Combining these result 
yields Eq. (4.20a) with coefficients appropriate for the 
case A >> qvp. Corrections are of 0((qvp/ A) 2 ). 



2. The functions /n, /12, and gu 

In order to calculate ipn, Eq. (4.18g), it is convenient 
to consider consider (pn(k—q, iSl n ), which can be written 
in a form similar to (p^(k,iCl n ), namely, 



where 



ipu(k - g,iQ. n ) = ip[i(k - q,iQ„) + <Pn ( fe _ Q,^n), 



(2), 



(k- q ,in n) = ^ w+w -- g(p) + A2 r ^ 



)-/(# + «,+ ) ; /(£+„,_)_/(£_„,+) 



+ 



iQ n — k ■ p/m c — w + — w- iQ n — k ■ p/m c + w + + w 



(C6a) 



(C6b) 



V$(k-q,iiln) = 



w+w- + g(p) - A 2 [ /(£ -«>-)- /(# - w+) + m + w-) - m q + + w+) 



W + W- 



iQ n — k ■ p/m c + w + 



W- 



iQ, n — k ■ p/m c — w + + w- 



. (C6c) 



r 



To evaluate these integrals we again assume A » qvp. In 
this limit 



(u> + w_ — g(p) + \ 2 )/w + W- 



(C7) 



and the integral again reduces to the ferromagnetic case. 
The result is 



<Pu(k- q,i£l n ) = <pu(q,i0) 



+ - 



_i / A y 

. 3 V2fc F / 



Shifting k by q, we find 

<p$(k,iSl n ) = ^(O.iO) 
N F \ 2 k q 1 



+ - 



3 (2/c F ) 2 3 V2fc F 



~2A~ 



~2A~ 



(C8a) 



. (C8b) 



For we need 

(w+w_ + g(p) - \ 2 )/w+w_ - (q • p) 2 /2m 2 A 2 . (C9) 

In contrast to Eq. (C4), this is quadratic in the inte- 
gration variable p. therefore is a stress correlation 
function. It has a hydrodynamic contribution of the same 
functional form as the Lindhard function, but only the 
transverse (with respect to fe) components contribute to 
it. The non-hydrodynamic contributions are again sub- 
leading compared to <fii ■ One finds 



<Pu ( k - <?> i^n) = <fu (k - q, i0) 



- ,J 2 h 



32 A Ifcl 3 A 



(ClOa) 
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and after shifting the momentum we obtain ten 

(2) /, -o ^ (2)// -n\ , at n Vf 1 fc i 
<p\{ (fe, ifi„) = <p K u ' ( k > l °)+ N F^ — -^ — ■ 

(ClOb) 

Using these results in Eqs. (4.18d, 4.18c) and (4.29c) we 

obtain Eqs. (4.20b, 4.20c) and (4.30a), respectively. ^(fc,^) = r}$(k,iil n ) + r}$(k,i{l n ), (Clla) 



The functions h. 



Finally, the function ij^i which is defined by Eq. (4.29e) 
and determines h^i according to Eq. (4.29d), can be writ- where 



p 



A x - w + + w_ 



W + W- 



w + — W- 



AV * — ' w + W- 
p 



m - w-} m + w+) f(e + w~) m ™+) 

i£l n — p ■ k/m c — w+ — W- 
/(£-«,_)-/(£-«,+) 



ifl n — p • k/ m c + w + + W- 



ifl n — p ■ k/m c + w + — w_ 
I 



in n — p ■ k/m c + w + — w_ 



(Cllb) 
(Cllc) 



Using the same techniques as above we find, in the limit r/^ is proportional to iVp(if2„/ \)(q-k) 2 jk\k 2 , and hence 

A ^> qvp, j g sm£ jj compared to rf^l by a factor of q 2 These 

Wfu -n > at i0 " rnn^ results yield Eq. (4.30b). 

%{{k,i9, n ) = -N F — . (CI Id) 
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